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ABSTRACT 

The correlation function £ (r) of matter in the non-linear regime is assumed to be deter- 
mined by the density profiles p(r) and the mass distribution n(M) of virialized halos. 
The Press-Schechter approach is used to compute n(M), and the stable clustering hy- 
pothesis is used to determine the density profiles of these Press-Schechter halos. Thus, 
the shape and amplitude of £(r) on small scales is related to the initial power spectrum 
of density fluctuations. 

The case of clustering from scale-free initial conditions is treated in detail. If n is 
the slope of the initial power spectrum of density fluctuations, then stable clustering 
requires that £(r) oc r~ 7 , where 7 is a known function of n. If halo halo correlations 
can be neglected, then p(r) oc r~ e , where e = (7+3)/2 = 3(4+n)/(5 + n). For all values 
of n of current interest, this slope is steeper than the value 3(3 + n)/(4 + n) that was 
obtained by Hoffman & Shaham in their treatment of the shapes of the outer regions 
of collapsed halos. Our main result is a prediction for the amplitude of the non- linear 
correlation function. The predicted amplitude and its dependence on n are in good 
quantitative agreement with iV-body simulations of self-similar clustering. 

If stable clustering is a good approximation only inside the half-mass radii of Press- 
Schechter halos, then the density contrast required for the onset of stable clustering can 
be estimated. This density contrast is in the range ~ 300 — 600 and increases with the 
initial slope n, in agreement with estimates from iV-body simulations. 
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1 INTRODUCTION 

Large-scale structure in the universe is thought to arise 
through the gravitational clustering of matter. In the nonlin- 
ear regime most of the dark matter is in bound halos which 
have separated out from the expanding background. There- 
fore the auto-correlation function, £(r), of the dark matter 
is closely related to the shapes of halos. In this paper we 
explore the consequences of this connection for the shape 
and amplitude of £(r) on small scales. 

There are two approximations that are commonly used 
to describe the evolution of gravitational clustering in the 
non- linear regime. The first is the assumption that, in this 
regime, the clustering is statistically stable. By this, one 
usually means that non- linear, virialized objects no longer 
participate in the expansion of the background Universe; 
they maintain their shapes in physical coordinates, so they 
shrink in comoving coordinates (e.g. Peebles 1965, 1980). 
The second approximation was first formulated by Press & 



Schechter (1974). They assumed that, on average, structures 
collapse spherically and then virialize, with initially more 
dense regions collapsing first, and less dense regions collaps- 
ing later. Further, they assumed that when they virialize, 
all halos have the same density (about 200 times the back- 
ground density at the time of virialization) whatever their 
mass. By applying these assumptions to an initially Gaus- 
sian density field, they derived the distribution of halos as a 
function of mass. In the Press-Schechter approach, the clus- 
tering is hierarchical, and virialized halos at a given epoch 
are progenitors of the halos that virialize at a later epoch. 
It is not obvious that these two idealizations, stable cluster- 
ing, and virialization at a fixed multiple of the background 
density at the time of virialization, are compatible. 

For example, consider a Press-Schechter halo with uni- 
form density (it has a 'tophat' density profile). Assume 
that at the time it virializes, say ti, it has mass Mi and 
density 178/9b(ti). This sets its radius Ri. In the Press- 
Schechter approach, at some later time t-z > ti, it will 
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have merged with other halos into an object of greater 
mass M2 > Mi. This more massive object will have density 
178 Pb{t2) < 178pb(ii)- On the other hand, if stable clus- 
tering is correct, then the core Mi within R\ will still have 
density 178 Pb(ii) so that the shell of mass (M2 — M1) around 
the core must be less dense than the core itself. Thus, if it 
had a tophat density profile at time ti, then the halo will not 
have a tophat density profile at the epoch ti - In other words, 
if Press-Schechter halos are to evolve consistently with the 
stable clustering assumption, then they cannot have tophat 
density profiles. This toy example illustrates that by requir- 
ing consistency between Press-Schechter and the stability 
approximation we can constrain the allowed density profiles 
of halos. 

In particular, it is known that the stable clustering ap- 
proximation allows one to constrain the shape of £(r) in the 
non-linear regime (e.g., Peebles 1974a; section 26 in Peebles 
1980; section 5.4 in Padmanabhan 1993). Section 2 uses this 
stable clustering shape of £(r) to constrain the density pro- 
files of Press-Schechter halos. The formalism of McClelland 
and Silk (1977) is then used to compute the amplitude of 
£(r) in the nonlinear regime. The predicted amplitude is 
compared with the results of N-body simulations. The den- 
sity profile inferred in Section 2 is different from the value 
that is obtained in other descriptions of spherical collapse, 
such as the Hoffman-Shaham model. Section 3 discusses the 
reasons for the difference. Section 4 studies the density con- 
trast required for the onset of stable clustering. It argues 
that only the cores of Press-Schechter halos are likely to 
have evolved consistently with both Press-Schechter and 
stable clustering. This has the consequence that the den- 
sity required for the onset of stable clustering should be an 
increasing function of n, where n is the slope of the power 
spectrum of initial fluctuations. Section 5 discusses the case 
of non-power law density profiles and summarizes the re- 
sults. 



2 THE CORRELATION FUNCTION IN THE 
NON-LINEAR REGIME 

In the highly non-linear regime, stable clustering may be a 
good approximation. This section combines the stable clus- 
tering hypothesis with the Press-Schechter model to com- 
pute the shape and amplitude of the correlation function 
£(r). We consider first some simple and general properties 
of£. 

2.1 Stable clustering and the density profile 

Consider the auto-correlation function of matter £(r) within 
a single halo. If the density profile of the halo is a power law, 



where 7 = 2e — 3 . 



(•3) 



p{r) = Ar % 
then £(r) is 

£(r) oc / d 3 s p(s) p(f + s) oc / dss 2 s _E |r + s\ 



(1) 



(2) 



which is just a convolution of the density profile with itself. 
When the initial profile satisfies 3/2 < e < 3, then the inte- 
gral converges to give a power law form for the correlation 
function (Peebles 1974b; McClelland & Silk 1977): 



If the density profile is a power law with slope e on small 
scales, but has a cutoff on large scales, then the correlation 
function is also a power law, with slope 7 given by equa- 
tion (^) on small scales, but with a cutoff on larger scales. 

Assuming that all matter is in halos of various masses, 
£(r) depends on correlations within halos (essentially the 
density profile), correlations between halos, and the distribu- 
tion of halo masses and radii. McClelland & Silk (1977) show 
that if all halos have the same mass and radius, and these 
halos are randomly placed in space (so that there are no 
halo-halo correlations), then the above result (equation ^|) 
still relates the shape of £ to the density profile p. They also 
show that, for randomly placed halos having a range of dif- 
ferent masses, but the same density profile, an integral over 
mass must be included. This integral affects the amplitude, 
but not the shape of £. (If the density profile has a mass 
dependent cutoff, then the position of the 'knee' in £ will 
depend on the number density of halos as a function of halo 
mass.) 

Of course, if the halos are not randomly distributed, 
then their correlations will change the shape of £. Since ha- 
los are likely to be correlated with each other, the relation 
(equation ^) between the shapes of £ and p would seem to 
be of limited applicability. However, if we restrict attention 
to small scales (the highly non- linear regime), then, for suf- 
ficiently small separations r, both members of each pair of 
particles will almost certainly be drawn from the same halo. 
This means that, for small separations, the shape of the cor- 
relation function is not affected by halo-halo correlations. If 
the halo-halo correlation function can be neglected, then 
we can use the formalism of McClelland & Silk (1977) to 
compute the shape and amplitude of £(r) in the non-linear 
regime. To do so, however, we must know the density pro- 
files of halos, and the number density of halos as a function 
of halo mass. 

In the highly non-linear regime, stable clustering is a 
good approximation, and it can be used to constrain the 
shapes of virialized halos. For clustering from scale-free 
initial conditions, Davis & Peebles (1977) show that the 
BBGKY hierarchy of equations admit self-similar solutions 
(Peebles 1980; Section 73). They show that, in the stable 
clustering regime, £(r) must be a power law: £(r) oc r~ 7 
with 7 = 3(3 + n)/(5 + n), where n is the slope of the ini- 
tial spectrum P(k) oc k n . This form for £(r), combined with 
equation (fa), yields a value for e in the stable clustering 



regime: 



7 + 3 _ 3(4 + n) 
2 (5 + n) '• 



for n > —3 . 



(4) 



Note that whereas e S c may be thought of as the slope of 
the two-point correlation function, subject to the constraint 
that one member of each pair of particles is certainly at the 
center of a halo, 7 is the slope of the correlation function 
when the locations of the members of each pair are uncon- 
strained. For < 7 < 3, equation (^J) shows that e > 7, 
which means that the correlation function centered on peaks 
is steeper than the unconstrained correlation function. This 
behaviour is similar to that of peaks associated with Gaus- 
sian random fields; the density profile of a typical spherically 
averaged Gaussian peak is steeper than the unconstrained 
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correlation function of the Gaussian field (Bardeen et al. 
1986). 

In the remainder of this paper we will use this density 
profile p(r) oc r -Eac to describe the shapes of virialized halos. 
The differences between this profile and that suggested by 
Hoffman & Shaham (1985) are discussed in Section 3. The 
effects of non-power law density profiles such as that pro- 
posed by Navarro, Frenk & White (1995) will be considered 
in Section 5. 



2.2 The correlation function 

The previous subsection argued that the McClelland & Silk 
formalism can be used to compute to shape and amplitude 
of the correlation function in the non-linear regime. To do 
so requires knowledge of the density profiles of halos, as well 
as the number density of halos as a function of mass. It 
argued that the stable clustering hypothesis could be used 
to determine the density profiles of the halos. This section 
uses the Press & Schechter (1974) model to determine the 
mass function. Then, it combines the stable clustering den- 
sity profile with the Press-Schechter mass function to com- 
pute the amplitude of the matter correlation function £(r) 
in the non-linear regime. 

In the stable clustering regime, we can use the McClel- 
land & Silk formalism to express £(r) as 



pjj M 2 X a (r)p(M,a)dMda 
[fj,jf Mp(M,a)dMda] 2 



(5) 



where p, is the total number density of halos, p(M, a) dM da 
is the probability that a halo has mass in the range dM 
about M and radius in the range da about a, and 



A„(r) 



u a {s) u a (s + r) ds, 



(6) 



where u a describes the average density profile of a halo 
of radius a (note the similarity to equation ^|). Here, 
M u a (r) dr — M, where M is the total mass of the halo. 
Thus, in this regime, £(r) is determined solely by the distri- 
bution of masses, radii, shapes, and density profiles of halos. 
Below, we will combine the Press-Schechter model with the 
stable clustering hypothesis to specify these distributions. 

In the Press-Schechter description of non-linear cluster- 
ing, all matter is assumed to be contained within spherical 
virialized halos, and the distribution of halo masses is de- 
termined by the initial conditions. For a scale free initial 
spectrum P(k) oc k n and fl = 1, the number density of 
halos with mass in the range dM about M is 



f n + 3 \ dM 

V 3 ) M 2 { ' 



(e.g. Efstathiou et al. 1988). The second scale, denoted r*, 
is the virial radius of an M* halo, and is defined below. 

The average density within each virialized halo is also 
known, so that the mass and radius of a virialized halo are 
related. For the halos described by equation ([?]) above, M = 
(47ra 3 /3) A n ip, where A n i = 178, and p is the background 
density as before. Thus, in the Press-Schechter model, the 
distribution of halo masses, radii, and shapes are specified, 
whereas the density profiles of these halos are not. So, in 
the Press-Schechter approach, the correlation function in 
the non- linear regime is 



/ M 2 X M (r)n(M) dM 
[jMn(M)dM] 2 



(9) 



where n(M) dM is given by equation (]7j). Notice that equa- 
tion (Q) implies that the denominator in equation Q is just 
p 2 . Moreover, notice that since A(r) (equation |^) is an ex- 
pression like equation (^) , equation Q is exactly of the form 
discussed in the previous subsection. That is, it is a sum 
over all the different masses of terms like equation (^). This 
means that the shapes of the correlation function and the 
density profile are related to each other (equation ^). 

To proceed, we need expressions for Aa/(?~), which 
in turn depend on the density profiles u a of the Press- 
Schechter halos. If we assume that these density profiles are 
power laws in distance from the halo center, then stable clus- 
tering requires that the density profile must be u a oc r~ e 
with e given by equation (|j|). Thus, the Press-Schechter 
model specifies the number density of halos as a function 
of mass, as well as providing a relation between the mass 
and radius of each halo, and, in the non-linear regime, the 
stable clustering hypothesis specifies the density profiles 
of the Press-Schechter halos. So, by combining the Press- 
Schechter model with the stable clustering hypothesis, we 
are able to compute the shape and amplitude of £(r). 

If we use the stable clustering value, e = 3(4+n)/(5+n), 
for the slope of the density profile of Press-Schechter halos, 
then we can compute \m{t). For —2 < n < 1, we know that 
e > 2, so 



Am(t-) 



47TO- 



= (— v 

V47ra 3 / r(e-2) 



\s — r\ 2 E — s 2 e l s 1 E ds, (10) 



/2 



for all < r < 2a, and Ajvf(r) = otherwise. Similarly, 
when n = —2, so e = 2, then 

ds 



A w (r) 



47T<T 2 r 



hi 



r/2 



— if < r < 2a, 

s 



(11) 



(Press & Schechter 1974; Lacey & Cole 1993). Here 

M» oc a 6/(n+3) (8) 

is a characteristic mass which grows as the universe expands. 
This characteristic mass defines two scales. The first scale, 
which we will denote tq, is defined as follows. If the variance 
in the mass contained in randomly placed cells of radius ro 
is 1.68 2 /2, then the mean mass contained in these spheres is 

M* = (47rro/3) p, where p oc a -3 is the background density If we absorb all constants into one term, then we can write 



and Am(0 = otherwise. 

In the limit of vanishingly small separations, 



A M (r) 
A M (r) 
Am(0 



Ajw(r) 



-(-) ■ 

ia z \rj 

I?)"' 



4:ira 3 
1.22 

47T<7 3 

1.58 

47T<7 3 



tt 2 /4 fa 
Aira :i 



5/2, 


n 


= 1, 


(12) 


2.4, 


n 


= 0, 


(13) 


9/4, 


n 


= -1, 


(14) 


2. 


n 


= -2. 


(15) 
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Ajw(r) = A n cr 7 3 /r 7 , so that the correlation function is 




where we have used the fact that 47rA n ip<r 3 /3 = M . (No- 
tice that this means that r 3 = 3M*/47rA n ip, so that r» is 
quite a lot smaller than the scale ro on which the variance is 
1.68 2 /2.) This, then, is the correlation function in the highly 
non-linear, stable clustering regime. 

It is interesting to compare this result with that which 
we would have obtained had all halos had the same mass, 
say Mhaio = f 3 M* . In this case, equation (^) shows that the 
correlation function is just £(r) = X n (47rA n i/3) (/ r*/r) 7 , 
for r < rhaio, so that the amplitudes are comparable if / 7 = 
r(ll+n/10+2n)/r(l/2). Note that the amplitude increases 
as / increases, and is comparable to the Press-Schechter 
amplitude when Mhaio is considerably less than M* (0.42 
and 0.125 times M», for n = 1 and n — —2, respectively). 



2.3 Comparison with simulations 

To illustrate the accuracy of the amplitude predicted by 
equation (^]), we can compare it directly with that mea- 
sured in the TV-body simulations of clustering from scale 
free initial conditions. A simple and reliable way to make 
the comparison is to use recently proposed fitting formu- 
lae for the non-linear £ which have been calibrated using 
high resolution TV-body simulations. These are based on the 
ansatz of Hamilton et al. (1991) which provides a universal 
relation between the linear and non-linear correlation func- 
tion (also see Gott & Rees 1975). This ansatz has recently 
been extended and refined (Nityananda & Padmanabhan 
1994; Peacock & Dodds 1994; Padmanabhan et al. 1995; 
Padmanabhan 1995; Jain, Mo & White 1995 - JMW). We 
shall use the fitting formulae proposed by JMW which take 
into account a simple dependence on the shape of the initial 
spectrum that is exhibited by the TV-body data. 

In the asymptotic non-linear regime the formulae are 
consistent with stable clustering and take the form (see 
equations (5a) and (6a) of JMW), 

£ NL (a,z)«f (i±^)- a4 C - L(a ,/)3/2 (17) 
where, 

I 3 « x 3 £ NL (a, x) ; t(x) = 4 / dy y 2 £(») . (18) 

x Jo 

Here, £(a,x) is the volume integral of the correlation func- 
tion over a sphere of comoving radius x = r/a (the sub- 
scripts refer to the non- linear and the linear regimes, re- 
spectively). Since J) oc Z _( - n+3 \ it is easy to see that 
the correlation function in equation ( |l6| ) has the same scal- 
ing as equation (p7|). 

The predicted amplitudes of the non-linear correla- 
tion function (equations [h] and ^) can be compared after 



Table 1. Comparison of the amplitude of the non-linear £ 



n 


N BC 




Fractional "Error" 





3.5 


3.1 


13 % 


-1 


5.1 


5.0 


2 % 


-2 


17.4 


13.7 


27 % 



some straightforward simplifications. For the initial spec- 
trum P(k) — Ak n , the expression in equation ( |l6| ) can be 
written as 

£ sc (a, r) = A TV sc (n) aT^J r" 7 , (19) 

where TV SC contains all the normalization factors. Also, equa- 
tion (^) simplifies to, 

6mw(a, r)=A TV jm w(n) r -7 . (20) 

So, the task of comparing the amplitude of £ to the TV- 
body calibrated formula simplifies to a comparison of the 
normalization constants TV sc (n) and TVj mw (n). The results for 
the cases n = 0, —1, and —2 are shown in Table 1. Note that, 
for the n = — 2 case, the agreement between our prediction 
and the TV-body data is better than indicated in the table 
because, as noted by JMW, their formula underestimates 
the non-linear amplitude of £ for this spectrum. 

We can also compare the amplitude predicted by equa- 
tion (^) directly with that measured in the TV body simu- 
lations of Efstathiou et al. (1988). Their Fig. 4 shows log 10 H 
where 3 = (x / X) 7 £(x / L) for —2 < n < 1. At late times, 
so that discreteness effects are less important, the ampli- 
tudes of their correlation functions scale as a 6 ^ 5+n ' ) . This 
scaling is consistent with that of equation (|lrj). At the fi- 
nal output time in their simulations log 10 H = —1.6,-1.4, 
and —0.9 for n = 1,0, and —1, respectively. If we set 
A n i = 178 and use the fact that, in their simulations, 
p = 32 3 /(aL) 3 and M» = C„ a 6/(3+n) , where C„ = 0.8, 0.71, 
and 0.53 (their eq. 15), then equation (jl6|) predicts that 
log 10 H = —1.6,-1.3, and —0.8, respectively. Thus, our pre- 
dicted amplitudes (equation |l(j) and those measured in the 
JMW and the Efstathiou et al. simulations agree to within 
about 20%. Some of this discrepancy is probably due to the 
fact that the Press-Schechter mass functions provide good, 
but not perfect, fits to the halo size distribution in the sim- 
ulations. We conclude that equation ( p^ ) provides a good 
description of the shape and amplitude of the correlation 
function in the highly non-linear regime. 



3 RELATION TO SECONDARY INFALL 
MODELS OF SPHERICAL COLLAPSE 

The calculation of the previous section assumed that Press- 
Schechter halos are spherical, with density profiles having 
slope e = 3(4 + n)/(5 + n). That is, the previous section 
provided a relation between the initial fluctuation spectrum 
and the density profiles of non-linear, spherical, virialized 
halos (equation^). It also showed that the correlation func- 
tion that corresponds to this choice for the relation between 
e and n is in good agreement with that measured in TV-body 
simulations. However, this relation between e and n can be 
tested more directly. 
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Recently, Crone et al. (1994) measured density profiles 
of virialized halos in their initially scale free simulations. 
Although they did not compare their measurements with 
the relation of equation (Q), their simulations show that it 
is remarkably accurate. Therefore, it is interesting to see if 
this relation can be derived directly from models of spherical 
collapse, rather than from the combination of the Press- 
Schechter model with the stable clustering hypothesis used 
in the previous section. 

It is well known that if the initial density profile of a 
halo is a power law in radius (a 'cone-hat', rather than a 
'tophat'), then the density profile of the collapsed halo is 
also a power law. If the initial density profile is p oc r~ a , 
then the collapsed halo has slope 



3a 
1 + a 



if a<2, 



(21) 



and e = -2 if a < 2 (Fillmore & Goldreich 1984). If one 
allows for non-radial orbits, then the restriction on a can 
be relaxed (White & Zaritsky 1992). This solution describes 
the shape of a single collapsed halo. If, instead, we start 
with a distribution of halos identified in an initially Gaussian 
random field, then the question is: What is the initial density 
profile of a typical halo? Once we know the answer, we can 
set this value equal to a in equation (j2^), and so calculate 
the typical density profile of collapsed halos. 

Hoffman & Shaham (1985) argue that the initial pro- 
genitor halos may be related to peaks in the initial density 
field. The density profiles of very high peaks have the same 
r dependence as the correlation function of the underlying 
field (Section VII in Bardeen et al. 1986). Since the corre- 
lation function of the initial (scale free) field is oc r~ f3+n ' , 
Hoffman & Shaham set a — 3 + n in equation ( [il] ) to ob- 
tain e = 3(3 + n)/(4 + n). This dependence on n is different 
from that in equation (tf). However, it is not clear that the 
Hoffman-Shaham value for a is the relevant one. This is be- 
cause, as the peak height decreases, the peak profile becomes 
steeper than the correlation function (Fig. 8 in Bardeen et 
al. 1986). As the universe expands, peaks of smaller and 
smaller height are able to collapse, so it is important to ac- 
count for this difference between the average peak profiles 
and the shape of the initial correlation function. 

The mean peak profile, averaged over all curvatures and 
orientations, depends, to a good approximation, on a sum 
that involves the initial correlation function and its second 
derivative (equation 7.10 in Bardeen et al. 1986). For ini- 
tially scale free fields, the correlation function is oc r _ ' 3+n ' 
so that its second derivative is oc r~' 5+n '. Thus, setting 
a — 4 + n should give some indication of the shapes of col- 
lapsed peaks, when the initial peaks are not arbitarily high. 
This value for a in equation (fell) gives e = 3(4 + n)/(5 + n) 
for — 2 < n < 1. While it is interesting that this relation 
between the non-linear density profile and the initial fluctu- 
ation field is the same as that obtained in the previous sec- 
tion using the stable clustering hypothesis, one would prefer 
a more rigorous derivation. 



4 THE ONSET OF STABLE CLUSTERING 

Section 2 used the fact that the shape of the correlation func- 
tion in the non-linear, stable clustering regime is related to 



the average density profile and mass distribution of virialized 
halos. It did not address the question of the scale on which 
stable clustering becomes a good approximation. The usual 
assumption is that stable clustering is accurate when the 
density contrast is on the order of that required for virial- 
ization: A n i ~ 200. On the other hand, iV-body simlulations 
show that the onset of stable clustering occurs at higher 
density contrasts as n, the slope of the power spectrum of 
initial fluctuations, increases (Jain 1995). This section uses 
the Press-Schechter description of clustering to show why 
this might be expected. 

Essentially, this section is motivated by the fact that not 
all Press-Schechter halos at the epoch ti evolve to become 
cores of halos by the epoch t2>t\. For example, one might 
reasonably expect that for a halo of mass M2 at the epoch t2, 
the most massive subhalo Mi < M2 at the epoch ti becomes 
the core. If the density profile of M2 at £2 has the slope 
obtained in the previous section, then some or all of the 
other subhalos that make up the mass (M2 — Mi) may have 
been destroyed as they merged to form the final halo. This 
disruption of subhalos as they merge to form larger halos has 
been seen in iV-body simulations (e.g. Fig. 2 in Efstathiou 
et al. 1988). These disrupted subhalos will not have evolved 
in accordance with the stable clustering hypothesis. So, only 
the (most massive?) progenitor subhalo in the core is likely 
to have evolved consistently with both Press-Schechter and 
stable clustering. This has the following consequence. 

It is usual to assume that stable clustering should be a 
good approximation on scales that are comparable to that 
of a typical Press-Schechter halo. This corresponds to scales 
on which the density is on the order of ~ 200 times the back- 
ground density. However, the discussion above suggests that 
the objects that are not in the core of a Press-Schechter halo 
may not satisfy the stable clustering hypothesis. Therefore, 
the relevant density contrast for the onset of stable clus- 
tering is that associated with the core of the halo, rather 
than of the halo itself. Suppose we assume that the core 
size is the scale associated with the typical size of the most 
massive progenitor of a given halo. Since the core is smaller 
(and more dense) than the halo itself, one might reasonably 
expect that stable clustering will only be a good approxima- 
tion once the density contrast is somewhat larger than the 
canonical value of ~ 200 or so. 

To be more quantitative, consider a halo with mass Mo 
at time to, and assume that a stable core was formed in 
its center when one of its progenitor subhalos had a mass 
of Mo/2. To estimate the density of this half-mass core, we 
must estimate the scale associated with this core, which, by 
the stable clustering assumption, is the same physical scale it 
had at the epoch at which it virialized. So, we must estimate 
the earliest epoch at which half the mass of the halo was first 
assembled into a virialized subhalo of mass Mo/2. There are 
two ways to do this. The first is to use the fact that 



P(Mi>M /2,ii|M ,io) =erfc 



V2 



(22) 



where P(Mi > M /2, ti |M , to) denotes the probability that 
a mass element of Alo was previously in an object of mass 
Mi > Mq/2 at the epoch ti, and where 



3(h) - S(t ) 



VK 2 (M,/2)-a ( 2(Mo))' 



(23) 
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where S(t) describes the (linear theory) overdensity required 
for collapse at time t (it decreases as the inverse of the expan- 
sion factor), and oq(M) is the variance of the linear density 
field when smoothed with a filter containing mass M (see 
Lacey & Cole 1993 for details). Then, setting P = 1/2 gives 
an estimate of 5(ti) = 8(tf) at the time tf of formation. 
Since erfc (x) = 1 — erf (x) = 1/2 when x fa 0.4769, this 
means that u>t = 0.4769 v2 = 0.67. This way of estimating 
the halo formation time has been used by Navarro, Frenk & 
White (1995). 

However, as Lacey & Cole (1993) note (and their Fig. 7 
shows), this method provides a biased estimate of the for- 
mation time. They show that a more precise estimate of the 
formation time is obtained by noting that, with this (half- 
mass assembled) definition of formation, different halos, each 
of mass Mo, may have formed at different times. So, one 
should compute the distribution of formation epochs, and 
then compute the mean, or the most probable value, of this 
distribution. For initially scale free fields, the most probable 
value of the formation epoch, tf, is given by 

! + «=(£) =1 + ^V^I(^) , (24) 

where /3 = (n + 3)/3, and Wt fa 0.75 (see Fig. 7 and 
equation 2.32 in Lacey & Cole 1993). This estimate differs 
slightly from the estimate of the previous paragraph, for 
which Wt fa 0.67. 

Equation (E4j) shows that the formation epoch is mass 
dependent. Fig! Ill shows this distribution of formation 
epochs for —2 < n < 1, with Wt = 0.7 (the curves for 
Wt = 0.67 and 0.75 are very similar to the one shown). Two 
features are obvious. Halos that are less massive than about 
10M* or so form earlier (at higher zt) for large n than for 
small n. On the other hand, halos that are more massive 
than this scale with n in the opposite sense. It is easy to 
see that both these features are a consequence of two facts. 
Namely, the distribution of Press-Schechter halo sizes be- 
comes broader as n decreases; when n = 1 most halos have 
the characteristic mass M*, but when n = —2, many halos 
have significantly different masses. So, if n = 1 initially, then 
halos that are significantly more massive than the character- 
istic mass at, say, the present time must have formed very 
recently (i.e., at low redshift), since they are extremely un- 
likely to have existed at early times when the characteristic 
mass was lower. In contrast, if n < —1 initially, then some 
of these more massive halos may easily have been formed 
at earlier times (higher redshifts). Furthermore, the growth 
rate of the characteristic mass is more rapid for small values 
of n (equation ^). Therefore, the asymptotic power law be- 
haviour shown in the figure (for small M/M* as a function 
of 1 + zt) is a direct consequence of equation (^). 

So, the way in which the epoch of formation of cores 
scales with n depends on the ratio M/M* . However, most of 
the mass is in halos with masses M < M, (the actual frac- 
tion is erf (1) = 0.84). This means that, for large n, most of 
the stable cores will virialize at earlier times than for lower 
n. By the Press-Schechter hypothesis, progenitor cores will 
have had an average density of ~ 200 times the background 
density when they virialized (at zt). This sets the scale of the 
halo. Stable clustering implies that, in physical coordinates, 
the size of the core does not change, so that, by to, its comov- 
ing density will have increased by the cube of the expansion 




M/M, 

Figure 1. The epoch of formation Zf for halos with mass M/M-*, 
for various choices of n. 

factor since the time tf when it virialized. So, on average, 
the average density of stable cores of Press-Schechter ha- 
los increases as n increases. Since stable clustering is only 
a good approximation in these cores, the density constrast 
on which stable clustering becomes a good approximation 
increases as n increases. 

We can estimate the amplitude of the correlation func- 
tion at which stable clustering becomes accurate as follows. 
Assume that, when averaged over all halos, the onset of 
stable clustering occurs at the half-mass radius of a halo 
with mass M = f s M* , where / was given at the end of 
section 2.2. Recall that this definition of / was chosen so 
that, if all halos had this same mass, rather than the Press- 
Schechter distribution of masses, then the r — > ampli- 
tude of the correlation function would be unchanged. The 
ratio of this half-mass radius to the virial radius of such a 
halo is 0.5 1//< - 3-e ' . On this scale, the correlation function is 
£ > \„ (47rA n i/3) 2 7/(3_E) which implies £ > 300 for n = -2 
and £ > 600 for n — 0. This increase of the requisite am- 
plitude of £ as n increases is in qualitative agreement with 
recent measurements of the onset of stable clustering in TV- 
body simulations for which £ > 200 for n — —2 and £ > 1000 
forn = (Jain 1995). 



5 DISCUSSION 

We have shown that requiring consistency between the sta- 
ble clustering hypothesis and the Press-Schechter approach 
suggests that the density profiles of Press-Schechter halos 
should be power laws with slope e = 3(4+n) / (5+n), where n 
is the slope of the initial fluctuation spectrum (equation ^). 
This value differs from that expected from secondary in- 
fall studies of the collapse of peaks associated with initially 
Gaussian fields (Hoffman & Shaham 1985). The reasons for 
this difference were discussed in Section 3. 

The Press-Schechter model, with the stable clustering 
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density profile, was used to compute the amplitude of the 
correlation function in the highly non-linear regime (Sec- 
tion 2). This amplitude was shown to be in good quantita- 
tive agreement with that measured in iV-body simulations. 
Moreover, the amplitude was shown to be a function of the 
slope n of the initial power spectrum. The n-dependence 
resembles that which was obtained recently by Jain et al. 
(1995) and Padmanabhan et al. (1995). 

Our calculation of £ (equation ||) rests on the assump- 
tion that, in this (small separation) regime, the correlation 
function is determined primarily by the density profiles of 
virialized halos. That is, the crucial assumption is that, for 
these small separations, most pairs are most likely to be 
drawn from the same halo, so that correlations between dif- 
ferent halos have a negligible effect on the shape of £. While 
this assumption makes intuitive sense, we have yet to show 
that it is accurate. We can do this as follows. 

Equation (g|) shows that the correlation function can be 
written as a sum over halos of different masses. To demon- 
strate that our approach is at least self-consistent, we must 
show that the correlation function is, indeed, determined 
primarily by halos having diameters that are larger than 
the separation scale r. So, consider the limit r — > of van- 
ishingly small separations. Then the amplitude of the corre- 
lation function is given by equation (161). Of course, in this 
limit, all halos are larger then the separation scale. However, 
as a crude approximation that should also give some indica- 
tion of the result on slightly larger scales, we will compute 
the fractional contribution to £ from 'small' halos, which we 
will define to be all those with M/M* < 0.2. Equation @ 
shows that this contribution is simply an incomplete Gamma 
function. Thus, the small mass contribution is 



7 (ll+n)/(10 + 2n),0.2 n+3/3 T (11 + n)/(10 + 2n) 



which is 11, 15, 20, and 24 per cent for n — 1,0,-1 and 
—2 respectively. For non-zero values of r, the corresponding 
integrals can be done numerically. They also show that £ is 
determined primarily by pairs within larger halos, so that 
our neglect of halo-halo correlations is justified. 

The other crucial ingredient in our calculation is knowl- 
edge of the density profiles of virialized halos. We assumed 
that these profiles were power laws, and determined the 
slopes of the power laws by requiring that they yield a corre- 
lation function with the shape required by stable clustering. 
However, recent work suggests that, while the profile shape 
may be described by an average power law slope that is 
consistent with equation (^), the density profiles of halos 
that form in iV-body simulations of gravitational clustering 
from scale free initial conditions are not simple power laws 
(Navarro, Frenk & White 1995; Cole & Lacey 1995; Tormen 
1995). Rather, the density profiles seem to be well fitted by 
a function of the form p oc l/[r(r + 6) 2 ], where the core ra- 
dius b depends both on the mass M and the shape of the 
initial power spectrum. 

We expect our formalism to work for these non-power- 
law profiles also because, except for the inner regions r ~ 
O.lrvir, where r v i r is the virial radius, they are reasonably 
well described by the power-law profiles we have been con- 
sidering. To illustrate this, Fig. ^| shows the non-power law 
profile p(s) oc l/[s(s + a) 2 ], where s = r/r v i r and the core 
radius is given by a — b/r v i r . The plot shows this profile for 




Figure 2. The density profile p(s) oc l/[s(s + a) 2 ], where s = 
f/fviri a = V r vin and b is the core radius, for a range of values 
of a. When s>0.1, the power law profile that corresponds to 
n = —2 (solid bold curve) describes the a = 0.25 profile quite 
well. The dashed bold curve shows the power law corresponding 
to the n = 1 case; when s > 0.1 it fits the a = 0.1 profile reasonably 
well. 



a range of values of the scaled core radius a. These values 
were chosen because iV-body simulations show that a ~ 0.25 
describes the n = —2 simulations quite well, whereas a is 
smaller (~ 0.15) when n = (Cole & Lacey 1995). For com- 
parison, the thicker curves show the corresponding stable 
clustering power law profiles for the two extreme cases stud- 
ied here, n — — 2 and n = 1. Notice that, in the outer regions 
where s>0.1, the n = —2 (solid bold) curve describes the 
a = 0.25 profile quite well, whereas the n = 1 (dashed bold) 
curve fits the a = 0.1 profile. 

Now, the correlation function in the non-linear regime 
is essentially an integral over the density profile (to give 
A m), followed by an integral over the mass function. For the 
choices of n and a shown in Fig. 1, the corresponding \m(t) 
curves for the power law and the non-power law cases are 
similar, on scales larger than about 0.1r v i r . This is hardly 
surprising, since Fig. |^ shows that the Navarro et al. and 
power law density profiles p{s) are similar in the outer re- 
gions where s>0.1 (which contain most of the mass of a 
halo), and A m is simply the density profile convolved with 
itself. 

Since, in our formalism, £ is simply the sum of many Am 
curves, this suggests that £ for these non-power law profiles 
should be similar to that for the power law profiles consid- 
ered in Section 2, at least for separations that are on the 
order of about 0.1 r* and larger. We have verified that this 
is indeed the case: for scales on which 1000, the pre- 

dicted amplitude of £ changes by less than about 10% for 
—2 < n < 1, if we use the Navarro et al. profile (with the 
values of a shown in Fig. 1) instead of power law profiles. 

On scales smaller than ~ 0.1 r*, it is necessary to take 
into account the dependence of the core radius on the mass 
of the halo. Simulations show that, typically, a decreases as 
M decreases. Navarro et al. (1995) argue that the qualitative 
form of this relation can be understood in terms of the for- 
mation times of Press-Schechter halos. Including this trend 
improves the agreement of our predicted £ with the result 
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from the Navarro et al profile. The asymptotic behavior of 
£(r) as r — > however, is more subtle, as it is very sensitive 
to the M — » behavior of the core radius a and the Press- 
Schechter n(M). Here we simply note that by requiring that 
£(r) have the stable clustering shape, we can obtain an in- 
dependent constraint on the relation between a, M, and n 
for these non-power law profiles. 

Since all of the results of this paper concern the stable 
clustering regime, Section 4 discussed the scale on which 
stable clustering should become a good approximation. It 
used the Lacey & Cole (1993) analysis of the merger histories 
of Press-Schechter halos to show that the onset of stable 
clustering occurs at higher density contrasts as the slope of 
the initial fluctuation power spectrum n increases. This is in 
qualitative agreement with iV-body simulations (Jain 1995). 

Our derivation of the amplitude of the correlation 
function in the stable clustering regime, using the Press- 
Schechter mass multiplicity function, can be extended to cos- 
mological models in which Q < 1. This is because the effect 
of Q, on the density profiles of halos has been calculated (e.g. 
Hoffman & Shaham 1985; Hoffman 1988; White & Zarit- 
sky 1992), as has the effect on the evolution of the Press- 
Schechter mass function (Lacey & Cole 1993) . Furthermore, 
recall that both the density profiles and the Press-Schechter 
multiplicity function depend strongly on the shape of the 
initial power spectrum. Thus, our calculation of £ shows 
explicitly how the non-linear correlation function depends 
on the shape and amplitude of the linear power spectrum. 
Therefore, if our approach is correct, then by requiring con- 
sistency between the shape of the mass multiplicity function 
and the shape of the non-linear correlation function, one can 
estimate the shape and amplitude of the initial perturbation 
spectrum. 
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